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1.  Introduction: 

During  the  contract  period,  the  principal  investigator  HJ.S.  Fernando,  graduate  student 
Michael  Levy,  undergraduate  student  Richard  Eijmbert  and  post-doctoral  associate  Sergei 
Voropayev  worked  on  several  problems  related  to  deep  convection  and  ice  formation  in  high 
latitude  oceans.  Several  new  phenomena  were  noted  and  useful  parameterizations  were  developed 
during  this  project  period.  Results  of  the  completed  studies  are  outlined  below.  Relevant 
publications  originated  from  the  P.I.’s  group  during  the  contract  period  are  listed  at  the  end  of  this 
report. 

Turbulent  Convection  from  Isolated  Sources 

An  experimental  study  was  carried  out  to  study  the  evolution  of  a  turbulent  buoyant  plume 
in  a  homogeneous  fluid.  This  study  differs  from  many  previous  investigations  carried  out  on  point 
plumes,  in  that  it  considers  finite  diameter  effects  and  the  nature  of  flow  near  the  source.  The 
results  have  applications  to  a  number  of  oceanic  and  atmospheric  convective  flows  wherein  typical 
horizontal  extent  of  forcing  D  is  much  larger  than  the  fluid  layer  depth.  Experiments  carried  out 
with  varying  source  diameters  D  and  source  buoyancy  fluxes  B0,  using  the  laser  induced 
fluorescence  and  particle  tracking  velocimetry  techniques  for  flow  diagnostics,  revealed  the 
following: 

(i)  The  initial  development  of  the  plume  during  the  time  t  <  1.2(d2  /  B0)1/3 occurs  as  if  the 
convection  is  horizontally  homogeneous;  i.e.,  the  depth  of  the  front  grows  as  h=0.35(B0t3)1/2,  the 
lateral  entrainment  flow  is  insignificant  and  the  frontal  propagation  occurs  by  the  engulfment  of 

ambient  fluid  by  turbulent  eddies.  The  entrainment  flow  could  be  detected  at  larger  times  t 
>1.2^D2  /  B0 )  ,  especially  after  the  plume  underwent  morphological  changes  and  achieved  a 

quasi-steady  state  at  a  time  te  =  1.8  (d2  /  Boy/3 . 

(ii)  The  plume  under  this  steady  state  consists  of  two  regions.  In  region  I,  0  <  z  <  zc,  the 
plume  width  shrinks  and  achieves  a  minimum  diameter  of  bm  «  0.55D  at  a  depth  of  zc  «  0.28D. 
In  region  II,  z  >  zc,  the  plume  expands,  much  the  same  way  as  a  point  plume,  and  at  larger  z  the 
source  diameter  effects  fade  away  thus  leaving  Bq  and  z  as  governing  parameters. 
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(iii)  Under  quasi-steady  conditions,  the  averaged  maximum  horizontal  velocity  in  region  I  was 
found  to  be 

Vx  -  0.66  (B0D)1/3, 

and  the  vertical  velocity  in  this  region  is  proposed  to  have  the  dependence 

Wc=2.7(B0D)1/3^£|. 

The  instantaneous  buoyancy  within  the  plume  shows  high  space-time  variability  but  stable 
statistical  averages  can  be  discerned.  For  z/D  <  0.1,  the  time-averaged  buoyancy  b(r,z) ,  spatially 
averaged  across  the  plume  width  at  a  given  z,  was  found  to  be 

b^lO/B^/D)173. 


(iv)  The  mean  centerline  velocity  and  width  averaged  buoyancy  measurements  in  region  II 
showed  the  following  for  larger  z/D  values: 


Wc  =1.2(B0z)1/3  for  z/D  >  0.3, 
and 


bav=3.7| 


(  p>2  ^ 
B0 


1/3 


for  z/D  >  0.4, 


indicating  the  absence  of  source  diameter  effects.  These  z  dependencies  have  similarities  to  that 
observed  in  turbulent  point  plumes. 


Experiments  on  the  Doming  of  Isopycnals  by  an  Upward  Suction  Velocity 

A  laboratory  study  was  carried  out  to  investigate  the  doming  of  the  density  interface  of  a 
two-layer  fluid  subjected  to  a  localized  upward  suction  velocity.  The  study  was  motivated  by  its 
possible  application  to  the  preconditioning  phase  of  deep  convection.  The  upward  flow  was 
induced  by  a  rotating  disk  placed  at  the  surface  of  the  fluid  system,  which  pumps  fluid  outward 
through  the  Ekman  layer  on  the  plate  and  induces  an  upward  velocity  near  the  center  of  the  plate  to 
compensate  for  the  radial  flow.  It  was  found  that  the  isopycnals  are  domed  when  the  Richardson 
number  Ri,  based  on  the  buoyancy  jump  across  the  interface  Ab,  the  upper-layer  depth  H,  and  the 
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suction  velocity  we,  Ri  =  AbH/  we2 ,  falls  below  a  critical  value  of  Ric=  1.6xl06.  As  the  doming 
increases,  the  density  interface  is  increasingly  exposed  to  the  turbulence  induced  by  the  disk  and  to 
the  secondary  flows  generated  at  the  tank  walls,  thus  causing  the  interface  to  recede. 

Development  of  a  Point  Plume  in  the  Presence  of 
Background  Rotation 

An  experimental  study  was  carried  out  to  investigate  the  evolution  of  a  point  plume  in  a 
rotating  fluid.  The  plume  delivers  a  specific  buoyancy  flux  B  at  the  surface  of  a  homogeneous 
fluid  layer  of  depth  H  subjected  to  a  background  rotation  rate  of  Q.  The  plume  descends  into  the 
fluid  layer,  deflects  at  the  bottom  surface  and  then  spreads  horizontally.  Depending  on  the 
experimental  conditions,  the  plume  can  be  affected  by  the  background  rotation  either  during  its 
descending  phase  or  the  horizontal  spreading  phase.  The  present  investigation  was  focused  on  the 
latter  case,  because  of  its  strong  geophysical  relevance.  The  major  findings  of  the  study  are: 

(1)  The  initial  descent  of  the  plume  occurs  as  if  the  background  rotation  is  absent,  and  the  time 
(t)  growth  laws  for  maximum  plume  width  bm  and  the  vertical  height  h  of  the  plume  front 
were  similar  to  that  of  non-rotating  fluids;  viz.,  bm  «  0.38(Bt3yM  and  h  =  1.7 ^Bt3^14 . 
After  traveling  a  distance  hci  =  3.3  (b/£23)1M  for  a  time  tci  =  2.4  Q-1,  a  measurable 
reduction  of  the  descent  rate  could  be  observed  because  of  the  impeding  influence  of 
rotation  on  the  lateral  entrainment  flow.  The  lateral  growth  of  the  plume  is  drastically 
reduced  at  the  time  of  tC2  ~  5.5  flf1,  possibly  because  of  the  influence  of  rotation  on  plume 
turbulence.  Thereafter,  the  plume  width  remains  at  an  approximately  constant  value  of  bc  = 
1 .4(B  /  Q3)174  but  the  vertical  descent  continues  until  the  plume  impinges  on  the  bottom  of 
the  tank. 

(2)  The  entrainment  flow  into  the  plume  is  modified  by  background  rotation  so  as  to  generate  a 
cyclonic  rim  current  in  the  upper  layers  surrounding  the  plume,  which  grows  in  thickness 
with  time  until  it  becomes  barotropically  unstable.  The  cyclonic  circulation  so  generated 
have  an  average  Rossby  number  of  um/2£2Dm  =  0.3-0.5,  based  on  the  maximum  velocity 
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um  and  the  corresponding  diameter  Dm .  For  H  <  hcl  case,  these  instabilities  are  observed 
only  after  the  plume  descends  to  the  bottom  of  the  tank  and  breaks  up  into  baroclinic 
eddies;  see  below. 

(3)  For  the  case  H  <  hcj,  the  plume  reaches  the  bottom,  deflects  and  spreads  laterally.  The 
lateral  spreading  of  the  deflected  plume  continues  in  an  approximately  axi-symmetric 
manner  until  the  radius  of  the  front  becomes  of  the  order  as  the  Rossby  deformation  radius, 
whence  the  front  is  decelerated  and  is  broken  up  to  form  anticyclonic  eddies.  As  such,  the 
observed  instability  appears  to  be  baroclinic  and,  as  stated  above,  preceded  the  instability  of 
the  upper-layer  cyclonic  flow.  The  wave  number  of  vortices  formed  depended  on  the  fluid 
depth  H,  in  addition  to  the  Rossby  radius,  signaling  a  strong  interaction  between  upper  and 
lower-layer  vortices. 

Although  in  our  experiments  the  flow  spread  over  the  bottom  surface,  no  stable  baroclinic 
vortex  was  formed  and  the  flow  always  disintegrated  into  several  vortices.  This  is  perhaps 
due  to  the  energetic  nature  of  the  deflected  flow,  which  possesses  additional  kinetic  energy 
acquired  during  the  plume  descent,  that  can  be  used  for  baroclinic  instabilities,  despite 
energy  dissipation  in  Ekman  layers. 

(4)  As  mentioned  above,  the  appearance  of  bottom  anticyclonic  vortices  precedes  instabilities 
of  the  upper  cyclonic  flow.  An  interesting  outcome  of  this  interaction  is  the  formation  of 
stable,  hetroclinic  vortex  pairs,  staggered  one  above  the  other.  The  interface  between  the 
layers  carrying  these  vortices  of  opposite  sign  domes  upward  and  the  vortex  pairs  translate 
slowly  in  the  background.  They  were  identified  as  having  close  similarities  to  the 
structures  known  as  “hetons.” 

(5)  Because  of  the  interest  of  applying  our  results  to  geophysical  convection  problems,  ranging 
from  small  hydrothermal  plumes  with  source  diameters  as  small  as  tens  of  centimeters  to 
deep  convective  processes  with  scales  of  the  order  tens  of  kilometers,  the  experiments  were 
extended  to  include  finite  diameter  ( d0)  effects  of  the  plume  source.  It  was  found  that  the 
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point-plume  approximation  describes  the  flow  fairly  well,  after  the  plume  descends  to 
depths  greater  than  about  12do-  • 

Some  Aspects  of  the  Decay  of  Convective  Turbulence 

The  decay  of  temperature  and  velocity  fluctuations  in  a  convective  turbulent  boundary  layer 
over  a  heated  surface  in  response  to  cooling  at  the  surface  was  studied  using  a  laboratory 
experiment.  The  time  scale  associated  with  cooling  was  greater  than  the  relaxation  time  scale  of 
turbulence,  whereupon  the  convective  layer  responds  to  the  changes  of  thermal  forcing  in  a  series 
of  quasi-stationary  steps.  It  was  found  that,  when  the  cooling  rate  \dTJdt\  was  maintained 

constant,  the  decay  times  of  turbulent  velocity  and  temperature  fields,  from  the  initiation  of 
cooling,  scale  with  AT0/\dTs /dt\,  where  A T0  is  the  temperature  difference  between  the  surface  and 

the  mixed  layer,  irrespective  of  the  Prandtl  number.  Mechanistic  explanations  were  provided  to 
describe  the  experimental  observations.  Geophysical  relevance  of  the  results  was  placed  in  the 
context  of  the  deep-ocean  convection  and  the  atmospheric  convective  boundary  layer. 

Turbulent  Convection  in  Rotating  and  Stratified  Fluids 

Turbulent  convection  induced  by  heating  the  bottom  boundary  of  a  wide,  linearly 
(temperature)  stratified,  rotating  fluid  layer  was  studied  using  a  series  of  laboratory  experiments.  It 
was  shown  that  the  growth  of  the  convective  mixed-layer  is  dynamically  affected  by  background 

rotation  (or  Coriolis  forces)  when  the  parameter  R  =  (h2/£23/q0)2/3  exceeds  a  critical  value,  R1; 

♦ 

which  is  approximately  100.  Here  h  is  the  depth  of  the  convective  layer,  Q,  is  the  rate  of  rotation 
and  q0  is  the  buoyancy  flux  at  the  bottom  boundary.  When  R  >  R2,  where  R2  is  approximately 
300,  the  buoyancy  gradient  in  the  mixed-layer  was  profoundly  affected  by  background  rotation. 
Conversely,  when  R  <  R2,  the  buoyancy  gradient  in  the  convective  layer  is  independent  of  the  rate 
of  rotation  and  approaches  that  of  convection  in  non-rotating  fluids.  When  R  >  R2,  the 
entrainment  velocity  was  found  to  be  dependent  on  the  buoyancy  frequency  of  the  overlying 
stratified  layer,  the  rate  of  rotation  and  the  conventional  (Deardorff)  convection  velocity.  A  simple 
theoretical  formulation  for  the  rate  of  entrainment  was  derived,  which  was  found  to  be  in  good 
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agreement  with  the  experimental  results.  The  experimental  results  also  indicate  that  entrainment  in 
this  case  is  dominated  by  non-penetrative  convection. 

The  Transition  from  Density-Driven  to 
Wave-Dominated  Isolated  Flow 

An  isolated  fluid  mass  travelling  horizontally  in  a  stratified  layer  is  a  phenomenon  that 
resembles  both  a  detached  gravity-current  head  and  a  strongly  nonlinear  solitary  wave.  Using 
laboratory  experiments,  we  examined  the  transition  in  time  form  a  regime  in  which  the  flow  is 
density  driven,  to  one  in  which  it  is  wave  dominated.  A  simple  means  of  creating  an  isolated  flow 
that  exhibits  both  density  and  wave  effects  was  achieved  by  dropping  a  thermal  into  a  linearly 
stratified  layer.  This  transitional  regime  is  called  an  ‘isolated  propagating  flow’.  Scaling  were 
developed  that  identify  parameters  for  which  the  transitional  regime  occurs.  Particle-tracking 
studies  revealed  the  vertical  flow  structure.  There  was  an  upper  zone  that  is  wave-  dynamical,  and 
a  lower  zone  in  which  transport  of  mass  occurred.  The  transported  mass  slowly  leaked  out,  until 
the  phenomenon  resembles  a  weakly  nonlinear  solitary  wave.  In  the  ocean,  cracks  in  the  ice  cap 
(polar  leads)  cause  similar  flows  impacting  the  thermocline. 

Ice  Formation  due  to  Surface  Freezing  of  a  Two-Layer  Fluid 

A  series  of  experiments  aimed  at  understanding  the  processes  associated  with  surface 
freezing  of  a  two-layer  fluid  was  carried  out  in  a  laboratory  tank.  The  flow  configuration  consists 
of  a  layer  of  cold,  salty  water  overlaying  a  relatively  deep  bottom  layer  of  warm,  saltier  water.  This 
situation  is  common  in  high-latitude  oceans  during  periods  of  rapid  ice  formation.  The 
experiments  were  conducted  in  a  tank  with  well-insulated  side  and  bottom  walls,  placed  in  a  walk- 
in  freezer  with  air  temperatures  from  -12  to  -20  °C.  A  system  of  thermocouples  was  used  to 
measure  the  temperatures  at  fixed  levels  in  water,  ice  and  air.  Microscale  conductivity  and 
temperature  probes  were  used  to  obtain  vertical  profiles  of  temperature  and  salinity  in  water.  In 
general,  when  external  fluxes  of  heat  and  salt  are  absent,  such  a  system  enhances  static  stability,  in 
the  sense  that  the  net  density  difference  between  the  layers  increases  with  time.  When  external 
fluxes  of  heat  (because  of  surface  cooling)  and  salt  (rejected  during  ice  formation)  are  applied, 
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however,  this  fluid  system  may  become  unstable  and  overturning  of  fluid  layers  is  possible.  In 
addition,  heat  transport  from  the  warmer  bottom  layer  to  the  colder  upper  layer  may  be  important, 
possibly  leading  to  the  reduction  in  the  rate  of  ice  formation  compared  to  that  of  a  homogeneous 
fluid  with  temperature  and  salinity  identical  to  the  upper  layer.  Descriptions  of  such  physical 
processes  were  advanced  using  laboratory  experiments,  and  quantitative  measurements  of  salient 
parameters  were  compared  with  the  predictions  of  a  theoretical  model  developed  to  explicate  the 
flow  evolution. 
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Velocity,"  submitted  for  publication.. 

Boubnov  B.M.  and  Fernando  H.J.S.,  "Experimental  Study  of  Convection  from  Isolated  Source 
of  Buoyancy,"  Submitted  (in  Russian)  to  Izvestia,  Fluid  Gas  Mechanics 

Boubnov  B.M.  and  Fernando  H.J.S.,  "Regimes  of  Convection  from  a  Local  Source  of  Buoyancy 
in  Rotating  Fluids,"  Submitted  (in  Russian)  to  Izvestia,  Atmospheric  and  Ocean  Physics. 

Boubnov  B.M.  and  Fernando  H.J.S.,  "Convection  from  a  Local  Source  of  Buoyancy," 
Submitted  (in  Russian)  to  Dokl.  RAS  (Proceedings  RAS) 

Conference  Proceedings 

Kit,  E.,  Fernando,  H.J.S.  and  Ching,  C.Y.,  Frequency  spectra  of  scalar  and  velocity  fluctuations 
at  entraining  stratified  interfaces,"  Sixth  European  Turbulence  Conference,  Lusanne;  In: 
Advances  in  Turbulence  VI,  Eds.  S.  Savrilakis,  L.  Machiels  and  P.A.  Monkewitz,  Kluwer 
Academic  Publishers,  595-597, 1996. 

Fernando,  H.J.S.,  "Discussion  on  the  Paper,  Conductivity  Measurements  in  the  Wake  of 
Submerged  Bodies  in  Density  Stratified  Media,"  Twenty-First  Symposium  on  Naval 
Hydrodynamics,  June  24-28,  Norwegian  University  of  Science  &  Technology,  Trondheim, 

1996. 

Fernando,  H.J.S.,  Peeters,  P.,  Robles,  A.  and  Niestadt,  B.  "Laboratory  Modeling  of 
Microbursts,"  1st  AIAA  Theoretical  Fluid  Mechanics  Meeting,  June  17-20,  New  Orleans, 
AIAA  Paper  96-2149,  1996. 

Kit,  E.,  Fernando,  H.J.S.  and  Strang,  E.  "Measurement  of  Zero-Mean-Shear  Turbulence  in 
Homogeneous  and  Two-Layer  Fluids,"  Second  International  Symposium  on  Hydrodynamics, 
Hong  Kong,  Dec.,  16-18,  In:  Hydrodynamics,  Theory  and  Applications  (Eds.  Chwang, 
A.T.,  Lee,  J.H.W.  and  Leung,  D.Y.C),  853-858,  1996. 

Fernando,  H.J.S.,  Colomer,  J.  and  Boubnov,  B.  "Laboratory  Studies  of  Convection  from 
Isolated  Sources,"  Proceedings,  Sixth  Asian  Congress  of  Fluid  Mechanics,  Madras,  857-859, 

1997. 

Strang,  E.,  Fernando,  H.J.S.  and  Kit,  E.  "Aspects  of  Turbulent  Mixing  at  Density  Interfaces," 
Proceedings,  Fifth  IMA  Conference  on  Stratified  Flows  -  "Mixing  and  Dispersion  in  Stably 
Stratified  Flows,"  Accepted,  Oxford  University  Press,  1997. 

Fernando,  H.J.S.  "Convection  in  Rotating  Fluids  and  Its  Geophysical  Relevance,"  Proceedings 
of  the  Konstantin  Federov  Memorial  Symposium,  St.  Petersburg.,  UNESCO  Volume  #  x,  in 
Press. 

Lozovatsky,  I.D.  and  Fernando,  H.J.S.  "Turbulence  in  Stratified  Patches  on  a  Shallow  Shelf," 
Proceedings  of  the  Konstantin  Federov  Memorial  Symposium,  St.  Petersburg.,  UNESCO 
Volume  #  x,  in  Press. 
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Conference  Presentations 


Ching,  C.Y.  and  Fernando,  H.J.S.  Modeling  of  Leads  as  Line  Buoyancy  Sources,  Am. 
Geophys.  Union,  Ocean  Sci.,  EOS,  76(3),  194,  1996 

Fernando,  H.J.S. ,  "Formation  of  Density  Staircases  in  Stratified  Fluids:  A  New  Mechanism," 
Bull.  Am.  Phys.  Soc.  41(9),  1820,  1996. 

Lozovatsky,  I.  D.,  and  H.  J.  S.  Fernando,  "Patchiness  and  Mixing  at  the  Shelf,"  The 
Oceanographic  Soc.  1996  Meeting,  Amsterdam,  1996. 

Fernando,  H.J.S.  "Rotationally  Dominated  Convective  Turbulence  and  Its  Relevance  to 
Convection  in  Oceans,"  29th  International  Liege  Colloquium  on  Ocean  Hydrodynamics. 

Fernando,  H.J.S.  "Modeling  of  Geophysical  Flows,"  Invited  Speaker,  Moscow'97  Conference 
on  Stability  and  Instabilities  of  Stratified  and  Rotating  Flows,"  Moscow,  1997. 

Fernando,  H.J.S.  "Laboratory  Experiments  on  Stratified  Shear  Flows,"  Keynote  Speaker, 
International  Workshop  on  Mixing  in  Geophysical  Flows,  Barcelona,  Spain,  1997. 

Voropayev,  S.I.  and  Fernando,  H.J.S.  "Laboratory  Experiments  on  the  Surface  Freezing  of 
Stratified  Salt  Water,"  Am.  Geophys.  Union,  Ocean  Sci.,  EOS,  79(1),  152,  1998 

Fernando,  H.J.S.,  Chen,  R-R.  and  Ayotte,  B.A.  "Plume  Dynamics  in  Rotating  Fluids,"  Am. 
Geophys.  Union,  Ocean  Sci.,  EOS,  79(1),  184,  1998 

Fernando,  H.J.S.  "Turbulence  and  Mixing  in  Environmental  Flows,"  Invited  Paper.  "Hopkins 
Conference  in  Environmental  Fluid  Mechanics,"  8-9,  April  2-4,  Baltimore,  MD  21218,  1998. 

Voropayev,  S.I  and  Fernando,  H.J.S.  "Laboratory  Experiments  on  the  Surface  Freezing  of 
Stratified  Salt  Water,"  30th  International  Liege  Colloquium  on  Ocean  Hydrodynamics,  May  4- 
8,  56. 

Fernando,  H.J.S.  "Convection  in  Rotating  Fluids  and  its  Relevance  to  Deep  Ocean  Convection," 
Federov  Memorial  Symposium  on  Oceanic  Fronts  and  Related  Phenomena,  May  18-22,  St. 
Petersburg,  Russia,  49,  1998. 

Fernando,  H.J.S.  "Stratified  Turbulence,"  Invited  Pedagogical  Lecture.  IUTAM.  IAMAS/IUGG. 
GTP/NCAR  Workshop  on  "Developments  in  Geophysical  Turbulence,"  Boulder,  CO  June 
16-19,  1998. 

Fernando,  H.J.S.  Entrainment  and  Mixing  in  Stably  Stratified  Flows,  Invited  Lecture. 
GTP/NCAR  conference  on  "Observations,  Experiments  and  LES"  Boulder,  CO  Aug.  12-15, 
1998. 

Fernando  &  Levy,  M.  "Convection  in  Rotating  Stratified  Fluids,"  Bull.  Am.  Phys.  Soc.  xx(ll), 
xxxx,  1998. 

Lucas,  D.L.,  Colomer,  J.,  Fernando,  H.J.S.  and  McCulley,  J.  "Urban  Heat  Islands:  Flow  and 
Temperature  Structure,"  Bull.  Am.  Phys.  Soc.  xx(ll),  xxxx,  1998. 

Mahalov,  A.,  Berman,  N.S.,  Fernando,  H.J.S.,  Yu,  F.  and  Pardyjak,  E.  "Stable  Layers  in  the 
Atmospheric  Boundary  Layer,"  Bull.  Am.  Phys.  Soc.  xx(ll),  xxxx,  1998. 
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Fernando,  HJ.S.  "Laboratory  Modeling  of  Deep  Convective  Processes,"  Labrador  Sea  Deep 
Convection  Meeting,  Big  Fork,  Montana,  Sep  9-11,  1988. 
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